Ignition and Front Propagation in Polymer Electrolyte Membrane Fuel Cells 
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Water produced in a Polymer Electrolyte Membrane (PEM) fuel cell enhances membrane proton 
conductivity; this positive feedback loop can lead to current ignition. Using a segmented anode fuel 
cell we study the effect of gas phase convection and membrane diffusion of water on the spatiotem- 
poral nonlinear dynamics -localized ignition and front propagation- in the cell. Co-current gas flow 
causes ignition at the cell outlet, and membrane diffusion causes the front to slowly propagate to 
the inlet; counter-current flow causes ignition in the interior of the cell, with the fronts subsequently 
spreading towards both inlets. These instabilities critically affect fuel cell performance. 

PACS numbers: 82.47.Gh,82.47.-a,47.54.-r 



Fuel cells constitute a reliable and environmentally 
friendly alternative energy source. Beyond steady state 
operation, understanding of their startup and transient 
dynamic behavior is crucial to their use in variable load 
automotive applications 0, 0] • These transient dynamics 
are intensely nonlinear due to a positive feedback loop: 
we recently demonstrated that water generated in a poly- 
mer electrolyte membrane (PEM) fuel cell increases pro- 
ton transport exponentially, which "ignites" the current 
0, 0] . There is a strong analogy to the positive feedback 
loop leading to nonlinear dynamics such as ignition and 
spatiotemporal front propagation in exothermic chemical 
reactions. There, a product of the reaction (heat) en- 
hances the reaction rate by raising the temperature ||. 
In autohumidified PEM fuel cells, the positive feedback 
between membrane water activity and proton conductiv- 
ity is known to cause steady state multiplicity 6]. 

Here we present spatiotemporal nonlinear dynamic phe- 
nomena (current ignition and current density front prop- 
agation along flow channels) in a two-dimensional fuel 
cell. PEM fuel cells typically have flow channels that dis- 
tribute the fuel (hydrogen) across the anode and oxidizer 
(oxygen) across the cathode (Fig. f^l 7]. Longitudinal 
water gradients in the membrane produce a sharp current 
density front that propagates in time along the channel 
8]. To design and control variable load, dynamic fuel 
cell operation, the mechanism underlying such dynamics 
must be understood. 

Experimental. Fig. [3 shows a photograph of a sim- 
plified fuel with a segmented anode, that permits current 
profile measurements. The current through each segment 
of the anode and the voltage drop across the external load 
resistor are recorded as a function of time. The cell oper- 
ates at atmospheric pressure; the flow rates of hydrogen 
at the anode and oxygen at the cathode are maintained 
with mass flow controllers. 

The fuel cell is placed between two temperature con- 
trolled aluminum blocks. There are six segments to the 
fuel cell, each with an active area of 0.5 cm 2 . It is pre- 
conditioned with a 20 load at 60 °C and flow rates of 
3.5 seem H2 at the anode and O2 at the cathode for 8-12 
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FIG. 1: (Color online) Schematic of a Hydrogen-Oxygen 
PEM fuel cell including internal and external equivalent cir- 
cuit elements. Hydrogen molecules dissociatively adsorb at 
the anode and are oxidized to protons. Electrons travel 
through an external load resistance to the cathode, while pro- 
tons diffuse through the PEM under an electrochemical gradi- 
ent to the cathode. Oxygen molecules adsorb at the cathode, 
are reduced, and react with the protons to produce water. 
The product water is absorbed into the PEM, or evaporates 
into the gas streams at the anode and cathode. 



hours. The current in each segment was < 1 mA after 
preconditioning (extinction) . 

The external load was then reduced to 0.25 £1 keeping 
temperature and flow rates fixed. The current in each 
segment remained at < 1 mA for more than 4 hours, at 
which time 100 /zL of water was injected into the an- 
ode feed stream. The current response for each segment 
after the water injection is shown in Fig. [3] The injec- 
tion of the water "ignited" the fuel cell current. Before 
water injection the resistance of the dry membrane for 
proton conduction was very high (> 10 ktt cm- 2 ), limit- 
ing the current and hence the water production. The in- 
jected water was absorbed into the membrane, decreasing 
the membrane resistance to ~ 1 fl cm -2 ; this increased 
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FIG. 2: (Color online) The segmented anode fuel cell. The 
anode electrode was broken into six individual elements sep- 
arated by Teflon spacers. Current through each element was 
measured independently. A membrane-electrode-assembly 
employing 2 ETEK electrodes with carbon supported Pt cat- 
alyst and a Nafion f 15 membrane was placed between anode 
and cathode. 



the currents to ~ 100 mA and made the fuel cell self- 
sustaining. 

To explore the location and dynamics of ignition in 
the flow channels, the fuel cell was extinguished as de- 
scribed above. The temperature was then reduced to 25 
°C, decreasing the water removal, and the load resistance 
reduced to 0.25 Q, increasing the water production. Cur- 
rent and voltage were recorded every 100 s for a period 
of 20 hours after these changes. For the measurements 
reported here, the fuel cell was placed with the flow chan- 
nels running vertically. Co-current measurements were 
made with both H2 and O2 flows going from top to bot- 
tom; for counter-current measurements, the oxygen flow 
went from top to bottom and the hydrogen flow from 
bottom to top. 

The ignition and subsequent current distributions 
along the flow channels as functions of time are shown 
in Figs. Ufa) and^Jb). An induction period of hours was 
required before any significant current was measured; the 
length of this induction period depended on how dry the 
fuel cell membrane was before startup (i.e. before reduc- 
ing the load resistance.) For co-current flow (Fig^Ja)) 
ignition first occurred at the outlet of the fuel cell; the 
current in anode element 6, near the outlet, rose from < 
1 mA to ~ 100 mA over a period of 5 min. A current ig- 
nition front then propagated from element 6 to element 
1 over a period of 15-20 minutes; as the current front 
propagated towards the entrance of the flow channels, 
the current at the exits of the flow channels dropped. 
With counter-current flow ignition first occurred at ele- 
ment 3, at the interior of the flow channel (Fig. Efb)). 
From the center, the ignition fronts "fanned" outwards, 
but the highest current always occurred in the center of 
the flow channel. The induction period was also longer 
with counter-current flow than with co-current flow. The 
long-term, total steady state steady state current was 
roughly double for counter-current (as opposed to co- 
current) flow. 

Modeling. We have been able to capture the basic 
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FIG. 3: Fuel cell ignition by water injection. The fuel cell was 
preconditioned so as to dry the membrane; the current in each 
segment of the anode was then < 1 mA. At t = 0, 100 fj,L of 
water were injected into the anode (hydrogen) feed. After the 
water was absorbed into the membrane, the fuel cell currents 
rose dramatically (ignited), and remained so. The numbers 
refer to the currents in the segments in order - 1 is at the 
anode inlet and 6 is at the anode outlet. 



physics of ignition dynamics in a simple lumped differ- 
ential, single anode fuel cell |(| ; to capture localized igni- 
tion and spatiotemporal front propagation the model is 
extended to a series of differential elements. The key fea- 
tures of the model are the water inventory in the polymer 
electrolyte, the transverse proton conductivity from the 
anode to the cathode and the longitudinal water trans- 
port through the membrane. Water sorbed into the poly- 
mer ionizes sulfonic acid groups, facilitating proton trans- 
port; water sorption by the PEM is limited by the total 
number of sulfonic acid groups in the membrane. Water 
is sorbed into, or desorbed from, the membrane depend- 
ing on the balance between water evaporating into the 
gas flow channels and water produced by the fuel cell. 
The multilayered membrane-electrode-assembly in PEM 
fuel cells results in complex intra- and interlayer trans- 
port processes. Since the flow channels are much longer 
than the membrane thickness (5 cm vs. 0.0125 cm), we 
make the simplifying assumption that water activity in 
the membrane is in local equilibrium with water activity 
in the gas flow channel above it, and the only gradients 
are longitudinal along the channel. 

The water balance in each differential element j of the 
membrane is given by eq. co-current gas flow (J = 1 
to 6; j = and j = 7 denote the anode inlet and outlet, 
respectively; the counter-current case easily follows); the 
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FIG. 4: (Color online) A comparison of the experimental 
and computed currents for co-current and counter-current 
flow of hydrogen and oxygen in a segmented anode PEM 
fuel cell. The color scale is for current through each an- 
ode segment in mA. (a) experimental co-current (b) experi- 
mental counter-current (c) computed co-current (d) computed 
counter-current. For the simulations, the flow rate for H2 and 
O2 is 3.5 and 6.5 mL min -1 respectively and k m = 3 x 10~ 6 
mol s _1 . 



inventory is balanced by water produced (1/2 the proton 
current), water convected in the gas flow, and longitudi- 
nal water diffusion (described by a lumped mass transfer 
coefficient between differential elements) . Eq. [5] is an 
empirical fit to the number of water molecules associated 
with each sulfonic acid group, A, as a function of water 
activity a w in a Nation 115 membrane 0. 



N< 



so 3 



dX(j) 



{V A 



Vn) — 
C1 RT 



a w (j) 



da w (j) 

[F A (j-l)+F G (j-l)]P w (j-l) 

k m [a(j + l) + a(j 



IT 

-[F A (j) + F c (j)}P w (j) 
-2a(i)] 



1) 



(1) 



AO") = 14.9 a w (j)- 44.7 a 2 w (j) + 70 a 3 w (j) 

-26.5 40')- 0-446 a 5 w (j) (2) 

We assume that the total gas pressure Pt is fixed, and 
the local water activity in the membrane is in equilibrium 
with the local water partial pressure P w , i.e. P w {j) = 
a w (j)P° and P H2 (j) = Po 2 (j) = P T - P w (j); P° is the 
water vapor pressure at this temperature. The molar flow 
rates change along the flow channel as water is formed; 
the molar flows are given by Fa {j ) — F A (j — 1) — ^0) / 4.F 
and Fc{j) — Fc{j — 1); the subscript A(C) corresponds 
to the anode(cathode), i(j) is the local current and T 
the Faraday constant. Lastly, we assume that the local 



potential Vfc between the anode and cathode is dictated 
by the thermodynamic driving force as in eq. |21 This 
neglects interfacial potential drops, which results in the 
predicted currents being about 20% larger than found in 
real fuel cells. 



VW(j) = 1.23-||ln 



Pl(j)PoAj) 



[V] (3) 



Based on the equivalent electrical circuit, the differen- 
tial elements are electrically connected in parallel to each 
other. The voltage across the external load resistance 
thus depends on the total current produced by all ele- 
ments; the local current is given by eq. 0] The local 
membrane resistance, RmU), depends of the local water 
content in the membrane. For a Nation 115 membrane 
employed in this fuel cell the membrane resistance as a 
function of water activity is given by eq. 9J. 
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In a single, differential PEM fuel cell ignition occurs when 
the initial water content in the membrane is sufficient 
for water production to exceed water removal by convec- 
tion. Here ignition will occur when (and where!) the 
water production exceeds water removal, i.e. when (and 
where) the right hand side of eq. becomes greater than 
zero. Water production depends on the load resistance 
and the membrane resistance. A dry membrane has a re- 
sistance of 500 kSl cm -2 , limiting the current density to 
a maximum of 2.4 /iA cm" 2 . According to eq. ^the feed 
flow rates would have to be < 0.1 mL hr" 1 at a fuel cell 
temperature of 60 °C for water production to be greater 
than water removal and ignite the fuel cell. Absorption 
of 10 /iL cm" 2 of water into the membrane reduces its 
resistance to ~ 10 fl cm" 2 and the maximum current 
density is 100 mA cm" 2 , sufficient for water production 
to exceed water removal, and the fuel cell to ignite as 
shown in Fig. [3] 

The key elements that account for ignition are (1) an 
exponential dependence of proton conductivity in the 
PEM with membrane water content and (2) the dynam- 
ics of water uptake into the PEM. The location of igni- 
tion and front propagation are consequences of (1) con- 
vection of water produced towards the cell downstream, 
where it can accumulate and (2) diffusion of water up- 
stream through the polymer membrane. Figs. Etc) and 
(d) show the simulated current profiles for co-current and 
counter-current flow with T = 47 °C, R L = 5fi. The 
model captures both the ignition and the front propaga- 
tion, though it overpredicts the experimentally observed 
currents. 

The model is only semi-quantitative because it ne- 
glects finite water mass transfer rates into the membrane 
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and from the membrane into the gas phases. It also 
neglects the effects of condensing liquid water, hinder- 
ing gas transport from the flow channels to the mem- 
brane/electrode interface. More detailed models that in- 
corporate these effects can give quantitative fits to ex- 
perimental results; yet the added complexity does not 
significantly enhance physical understanding. 

For co-current flow the water produced upstream is 
conducted towards the outlet, where it slowly accumu- 
lates in the membrane; when the water content increases 
to the point where the local membrane resistance be- 
comes comparable to the external load resistance the 
current starts increasing rapidly, hydrating the mem- 
brane and causing ignition at the outlet of the flow 
channel. Upon ignition the water activity in the mem- 
brane approaches unity. Water is then transported up- 
stream through diffusion in the membrane itself, causing 
upstream propagation of the ignition. The model did 
not capture the observed eventual decrease in the down- 
stream current, after the ignition propagated to the inlet 
of the fuel cell. This decrease is due to condensing water 
accumulating in the cathode, inhibiting oxygen transfer 
to the catalyst and reducing the current. 

In counter-current flow, water formed at the cathode 
is convected towards the anode inlet in the cathode flow 
channel. It is also transported across the membrane to 
the anode, where it is convected towards the anode out- 
let. Water accumulates fastest towards the middle of the 
flow channels, resulting in an interior ignition point. The 
ratio of the flow rates between the anode and cathode 
affects ignition point location; a relative increase of the 
anode flow rate shifts ignition towards the anode outlet. 
After ignition, water starts to accumulate locally in the 
membrane. The transport of water through the mem- 
brane from high concentration at the middle of the flow 
channel towards the outlets results in the "fanning out" 
of the ignition fronts. 

Liquid water was observed leaving the flow channels 
30-40 minutes after ignition. Ignition takes the fuel cell 
from very low water activity to water activity of unity: 
condensation of liquid water in the flow channels. Grav- 
ity plays a key role in how such liquid water moves 
through the flow channels; the operation changes dra- 
matically if gas flow in the channel is counter to gravity 
driven liquid water flow. When the fuel cell was vertical, 
gravity caused the liquid to drain, and permitted good 
access for the reactants from the flow channels to the 
electrode/electrolyte interface. When the fuel cell was 
horizontal the initial ignition phenomena were similar to 
those reported in Fig. [3J after ignition, however, large 
fluctuations in the local current density appear to corre- 
late with water droplets exiting the cell. In the horizontal 
orientation, liquid water condensing in the flow channels 
could partially block flow. The liquid drops were pushed 
along the flow channels by the flowing gas, but in an ir- 
regular fashion, that gave rise to large fluctuations in the 



local current density. Transport of liquid water in the 
flow channels is not accounted for at our level of mod- 
eling. The model can thus capture ignition and front 
propagation, but will break down at longer times, when 
liquid water floods the cathode gas diffusion layer (GDL) 
and starts entering the flow channels. 

We have demonstrated how the water exponentially in- 
creases proton conductivity in polymer electrolyte mem- 
branes leading to ignition of the current in PEM fuel 
cells. Water can be injected to a fuel cell to ignite 
the current, just like a match can be struck to ignite 
a flame. Transport of water laterally, coupled with the 
exponential increase in proton conductivity produces cur- 
rent fronts that propagate along flow channels, just like 
flame fronts. The positive feedback loop between wa- 
ter production and increased proton conductivity of the 
electrolyte membrane in PEM fuel cells is analogous to 
exothermic chemical reaction ignition and flame propaga- 
tion, yet "water fans the flame" in PEM fuel cells! Front 
propagation depends on flow configurations creating dif- 
ferent front ignition and propagation patterns. Under- 
standing the parametric dependence and time scales of 
these phenomena is a vital component of fuel cell design, 
non-steady state operation and control. 
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